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Nucleon-nucleon scattering in spin-triplet channels is analysed within an effective field theory 
where one-pion exchange is treated nonperturbatively. Justifying this requires the identification 
of an additional low-energy scale in the strength of that potential. Short-range interactions are 
organised according to the resulting power counting, in which the leading term is promoted to 
significantly lower order than in the usual perturbative counting. In each channel there is a critical 
momentum above which the waves probe the singular core of the tensor potential and the new 
counting is necessary. When the effects of one- and two-pion exchange have been removed using a 
distorted-wave Born approximation, the residual scattering in waves with L < 2 is well described 
by the first three terms in the new counting. In contrast, the scattering in waves with L > 3 is 
consistent with the perturbative counting, at least for energies up to 300 MeV. This pattern is in 
agreement with estimates of the critical momenta in these channels. 



I. INTRODUCTION 



Over the last fifteen years, considerable effort has gone into trying to analyse nuclear forces using the systematic 
tools of effective field theory (EFT).^ The starting point was Weinberg's original proposal that these forces could 
be described within the framework of chiral perturbation theory (ChPT). This approach organises the terms in the 
effective Lagrangian or Hamiltonian according to powers of low-energy scales they contain. These scales, generically 
denoted Q, include momenta and the pion mass. In this "Weinberg" power counting, the leading terms of the nucleon- 
nucleon potential are one-pion exchange (OPE) and an energy-independent contact interaction, both of which arc of 
order . 

Weinberg also noted the enhancement of the nonrelativistic two-nucleon propagator near threshold and proposed 
that the leading terms in the potential should be iterated to all orders in order to generate nonperturbative effects, 
such as the deuteron bound state. This approach, referred to here as the "Weinberg-van-Kolck" (WvK) scheme, has 
been widely applied by van Kolck and collaborators and by many others.^ However, even with this enhancement, 
nucleon-nucleon loop integrals are of order Q. Although this is lower than the order expected in a relativistic 
theory, it means that each iteraction of the leading potential in a scattering equation raises the order by one power 
of Q. Hence the resulting amplitude should still be perturbatively expandable in powers of the scales Q. 

In order to justify treating the leading terms in the potential nonperturbatively, a further IR enhancement is 
needed, to promote them to order Q^^ . Such a promotion is only possible within a consistent power counting if we 
can identify additional low-energy scales in the nucleon-nucleon system. In the case of S'-wave scattering, the large 
scattering lengths provide such scales, and these lead to an EFT in which the leading, energy-independent contact 
terms are treated nonperturbatively @, [1, 01- At low energies, where the finite-range of OPE is not resolved, the 
resulting expansion of the potential is simply the effective-range expansion Q . 

Although a similar systematic justification for the iteration of OPE was not provided, the WvK scheme has been 
successfully used to describe a variety of few-nucleon systems and their interactions. Nonetheless, its validity in the 
^51-^-01 channel has been questioned 10] and, more recently, several groups have observed that Weinberg power 
counting can break down for nucleon-nucleon scattering in spin-triplet channels with nonzero orbital angular momen- 
tum [ill [T^. [Tsj.^ In particular, Nogga, Timmermans and van Kolck [ll[ find that the leading contact interactions 
can be substantially promoted in channels where tensor OPE is attractive, although this conclusion does depend on 
the choice of cut-off, as stressed by Epelbaum and Meissner p^ . 

To establish a quantitative form for this new power counting, we need first to identify a low-energy scale that 
would justify iterating OPE, and then to analyse the scale dependence of the associated short-range interactions. 



^ For reviews, from various points of view, see Refs. [3, 

^ Examples of successful applications can be found the reviews cited in footnote 1. 

^ Closely related observations can be found in the work of Pavon Valderrama and Ruiz Arriola ,14,1 . 
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The renormalisation group (RG) [l^ provides the natural tool for such an analysis. In its Wilsonian version, it has 
been applied to two-body scattering by short-range forces, showing that the effective-range expansion is based on a 
nontrivial fixed point of the RG flow [l^ . Distorted- wave methods have been used to extend the approach to systems 
with known long-range forces [3 ■ 

Once a factor of 1 /Af „ has been divided out of the Hamiltonian, the strength of the OPE potential can be expressed 
in terms of the momentum scale 

2 

777 

= » M 290 MeV. (1) 

/it jvjv-'"n 

where /ttmn is the pseudovector ttN coupling constant. In the chiral limit, it can be written 

A. ^ (2) 

where Mjv the nucleon mass, is the axial coupling of the nucleon and -FV the pion decay constant. In strict 
ChPT, At is therefore a high-energy scale, built out of ^ttFtt and Mj^. None-the-less its numerical value is small, only 
about twice tti^. As a result, perturbative treatments of OPE (as advocated by Kaplan, Savage and Wise 7]) fail to 
converge or converge only slowly [l^, [l^ . This suggests that we should explore the consequences of identifying Att as 
a low-energy scale, counting it as of order Q. Since A^r is proportional to 1/Afjv, this can be thought of as a concrete 
version of Weinberg's suggestion that 1 /Afjv should be treated as if it were of order 

In Ref. [12j , I applied a renormalisation group analysis to the short-range potential in the presence of tensor OPE. 
This made use of the distorted waves (DW's) of a potential (the chiral limit of tensor OPE). In the resulting 
power counting, the leading short-range potential is of order Q~^^^, independently of the orbital angular momentum. 
This is quite different from Weinberg counting, where the leading term in the L-th partial wave is of order Q^^. 
Subleading terms containing powers of the energy appear at orders Q^^^, Q^^^, and so on. This promotion of short- 
range terms confirms the numerical observations of Nogga, Timmermans and van Kolck [Tlj and makes quantitative 
the new counting proposed there. The terms in the resulting potential can be directly related to a DW Born expansion, 
similar to that in Refs. pTlflst. 

The validity of this counting does depend on the energies considered, since in each channel there is critical momentum 
above which waves penetrate the centrifugal barrier and reach the region where the singularity dominates. The 
analyses of Refs. [ill, [13 show that a nonperturbative treatment of OPE, and hence the new counting, is needed 
in the 5, P and D waves for momenta of order ttIt^. In contrast, waves with L > 3 do not probe the singularity 
until momenta of ~ 2 GeV are reached. In these higher partial waves OPE can be treated as a perturbation and 
short-distance interactions can be organised according to the usual Weinberg power counting. 

The results of the analysis of Ref. |12| were purely formal, leading to the power counting that governs the importance 
of the terms in the expansion of the short-range potential. In the present paper, I explore its practical consequences 
by analysing nucleon-nucleon scattering in spin-triplet channels, with an extension of the method applied to singlet 
channels in Ref. [l^l- For simplicity I consider only the uncoupled waves: '^Po,i, ^D2, '^F^ and ^6*4. The extension 
of the method to coupled waves such as '^Si-^Di or ^P2-^f2 is very similar in principle, but is technically more 
complicated because of the matrix nature of the equations. 

The RG analysis relies on the forms of the DW's at small radii, where they tend to asymptotic forms that are 
independent of energy. These waves are obtained by solving the Schodinger equation, as described in Sec. II. At small 
enough radii they show nonperturbative behaviour controlled by the l/r^ singularity of the tensor potential. In the 
case of waves with L < 2, this region extends out to about 1 fm. For lab kinetic energies up to 300 MeV the waves 
reach their asymptotic forms only for radii less than about 0.6 fm, and there they are dominated by the 1/r^ potential. 
This nonperturbative behaviour is present in waves with L > 3, but only for only for radii less than about 0.2 fm. In 
the range 0.2-0.6 fm they have the normal power-law forms associated with the centrifugal barrier. This confirms the 
expectations in Refs. [Ill Il2l | that low partial waves need the new power counting for energies in this range, whereas 
the higher waves can still be described perturbatively using Weinberg counting. 

I then use DW methods to "deconstruct" empirical scattering amplitudes by removing the effects of known long- 
range-forces. The residual amplitude can then be interpreted directly in terms of an effective short-range potential. 
This technique can provide a better indication of how well the known forces are able to describe the scattering, 
compared to simply plotting phase shifts. Such plots can be misleading since they tend to hide small differences in 
peripheral waves at low energies, which is just where the long-range forces should dominate. If the resulting potential 
still shows strongly nonlinear energy dependence at low energies, then this implies that long-range forces are still 
making important contributions. Short-range forces lead to a smooth energy dependence that can be expanded as a 
power series. 

As in the similar treatment of singlet channels [l^], I take several Nijmegen PWA's or potentials [2l|, to give an 
indication of the uncertainties involved in these analyses of the data. In Sec. Ill, I use them to construct scattering 
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amplitudes between the DW's of the OPE potential and from these I extract short-range potentials in the uncoupled 
spin-triplet channels. The resulting potentials show rapid energy dependence at low energies, indicating that important 
long-range physics is still present. 

The most obvious long-range forces that need to be removed next are two-pion exchange (TPE) and relativistic 
corrections to OPE. These appear at orders and (in Weinberg counting). I use here the forms of the TPE 
potentials given in Refs. [H, [13] and the corresponding order-Q^ correction to OPE [23|- At this order, there is 
also a 77r-exchange potential, calculated in Ref . [25[ . These potentials can all be subtracted perturbatively using the 
DWBA. The residual short-distance interactions shown in Sec. IV are consistent with the new power counting in the 
■^Po.i and waves. In higher waves, ^Fj, and ^G\^ the uncertainties in the Nijmegen PWA's make it hard to draw 
very strong conclusions but the residual short-range potentials are smaller after removal of TPE, and similar to those 
in the singlet channels ^20.] . 



II. DISTORTED WAVES 



The radial Schrodinger equation that describes the relative motion of two nucleons interacting through the long- 
range OPE potential is 



1 



d^ 

dr2 
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where the central piece of the lowest-order potential is 

1 



and its tensor piece is 
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Here the tensor operator, 5*12 = 3(cri • f )((T2 • f) — (Ti ■ (T2, takes the value +2 in the uncoupled '^Pi, ^D2, ■ ■ ■ channels, 
and —4 in the '^Pq channel. The isospin factor, ti • X2, is -1-1 for channels with odd L and —3 for even L. The on-shell 
momentum in the centre-of-mass frame, denoted by p, is related to the lab kinetic energy, T, by T = 2p'^/M,^. 

At small enough radii, all the solutions in a given channel tend to a common, energy-independent form which is 
determined by the term of the tensor potential and the centrifugal barrier. This form can be found by 
solving Eq. ([3]) at zero energy in the chiral limit (m7r = 0) where it can be written 



dr^ r dr 
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Here I have specialised to the case of uncoupled triplet channels and I have introduced the length scale 
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The solutions in this limit can be expressed in terms of Bcssel functions of order 2L + 1. This can be seen by defining 
the variable x = \/\PLj\/r and the function (l){x) = x^^^'^^i:,{\(3lj\/x'^), so that the equation becomes 



dx2 
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[x) = 0, 



(8) 



where the plus sign applies to channels with even J (where the tensor potential is attractive) and the minus sign to 
odd J. 

The solutions in the attractive channels are oscillatory: 
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where J2L+1 and I2L+1 denote the regular and irregular Bessel functions 

-1/4 



In the limit r ^ 0, these solutions tend 
They depend on a free parameter a. This 



to the WKB form of a sinusoidal function of 2y/\PLj\/r times 
angle fixes the phase of the small-r oscillations or, equivalently, it specifies a self-adjoint extension of the original 
Hamiltonian (see Refs. [13. [Tsj for further references). This is necessary since both Bessel functions give acceptable 
solutions for small r. There is a redundancy between a and the leading term in the effective short-range potential 
since both have the effect of fixing the phase of the wave function for small r. In waves where the scattering is weak, 
it is simplest to set a = and use the potential to represent short-distance physics. This leads to solutions to Eq. ([6|) 
that grow like for large r. In channels where the OPE potential can be treated perturbatively, this allows the 
waves to match on to their usual short-distance forms at larger radii where the centrifugal barrier dominates over the 
tensor potential. 

The special choice a = 7r/2 gives a solution to Eq. ^ that decays like r^'^^+^' for large r. Imposing this boundary 
condition on the full OPE problem would lead to a wave that was very large inside the attractive well of the 
potential. This would correspond to a system with a low-energy bound state or resonance. Since none of the NN 
channels with L > has such a low-energy state, values of a close to tt/2 should be avoided. In practice, a = is a 
good choice for all but one of the waves studied here. 

In the repulsive channels, the solutions are given by modified Bessel functions, and the regular one has the form 
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This vanishes exponentially with 2y^\(3Lj\/r as r ^ but, like Eq. grows as for large r. 

It is convenient to normalise these solutions so that, as r increases, they match on to the expected short-distance 
behaviour of the free solutions, jL{p'r)/p^ ■ (Since the asymptotic solutions are defined at zero energy, I have divided 
out the energy dependence from the spherical Bessel functions.) With this normalisation, 
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the asymptotic solutions become 
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Fig-fflshows the solutions to the full Schrodinger equation ([3]) for the channels '^Pq, ^Pi, and "^04 at lab kinetic 
energies of 5 and 300 MeV, and compares them to the zero-energy, chiral- limit solutions of Eq. (|12p . To make this 
comparison easier, I have also divided out the energy dependence from the full solutions. The waves shown for the 
attractive channels all have the short-distance phase a = 0. 

Although there is still some energy dependence in their normalisations, the shapes of the solutions in these channels 
reach their energy-independent forms for radii smaller than about 0.8 fm. The waves in channels where the tensor 
potential is attractive all oscillate for small enough radii. In the lowest waves, '^Pq and '^D2, this behaviour covers the 
whole energy-independent region. Indeed the first node of the ^Pq wave function appears at the edge of or, depending 
on the choice of a, within the domain of the low-energy EFT. Moreover the normalisation of the short-distance wave 
functions shows significant energy-dependence, beyond the p^ expected from the free solutions. 

In contrast, the ^Gi wave becomes oscillatory only for radii smaller than about 0.25 fm, beyond the scope of 
our EFT. In the rest of its energy-independent region it has power-law behaviour as expected from the centrifugal 
potential. More importantly from the EFT point of view, the energy dependence of the short-distance normalisation 
is almost entirely given by the p^ factor expected for a free solution. These results are consistent with the estimates 
in Ref. based on the chiral limit of the tensor potential. There, the critical momentum scales for the breakdown of 
perturbation theory were found to be of the order of to^ or in the S, P and D waves, and much larger (~ 2 GeV) 
in the F waves and above. 

The and waves have low-energy bound states or narrow resonances if a is taken to be close to n/2. The 
phase shifts can depend significantly on the choice of a in this region which, in the ^D2 case, is roughly 7r/4 a < Sir/ A. 
The high centrifugal barrier in the ^G^ wave means its phase shift is very weakly dependent on a, outside a very 
narrow band around 7r/2. In contrast the ^Pq phase shift shows a strong dependence on a over the whole range 
< a < TT. 
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The wave functions in the repulsive channels '^Pi and '^^3 both have qualitatively similar forms, and so only one 
is shown in Fig. [TJ They switch from power-law behaviour at larger radii, where the centrifugal barrier dominates, 
to exponential at small radii, where the potential wins out. In both cases the energy dependence of the short- 
distance normalisation is quite well described by the free form, at least for energies below about 300 MeV. At 
higher energies than considered here, 500 MeV or above, the ^Pi normalisation does develop a much stronger energy 
dependence (whereas the '^^3 does not). These numerical observations indicate that the effect of the finite pion mass 
has been to shift the critical momentum scale in ■^Pi channel to a somewhat higher value than the estimate in Ref. T^. 
From the wave functions alone, it is thus not clear whether tensor OPE is better treated nonperturbatively in this 
channel. 



III. ONE-PION-EXCHANGE EFFECTIVE POTENTIAL 

The RG method developed in Refs. [l3| shows that the terms in the short-range effective potential are directly 
related to an expansion of the DW if-matrix in powers of the energy and other low-energy scales. A DW approach 
like this is obviously needed in channels where OPE must be treated nonperturbatively. I will also use it in the 
channels where a perturbative treatment of OPE would be valid, since it provides a convenient alternative to fourth- 
order perturbation theory. The scattering between DW's of the long-range potential can be described by a reactance 
matrix K. In the uncoupled channels, its on-shell matrix element is related to the difference between the observed 
and OPE phase shifts by 

K{p) - tan(^(5pwA(p) - (JopEb)) , (13) 

where <5pwa(p) is the empirical phase shift (from the Nijmegen group's PWA or one of their potentials that have been 
fitted directly to data) and (5ope(p) is that obtained from the solutions of Eq. ([3]). None of the channels examined 
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here contains a low-energy bound state or resonance, and so the residual scattering can be represented by an EFT 
expanded around a trivial fixed point. The amplitude K{p) is then given by the distorted-wave Born approximation 
(DWBA): the matrix element of the short-range interaction between the DW's of the long-range potential. 

Since these DW's are either vanishing or singular as r ^ 0, an ordinary 5-function at the origin cannot be used to 
represent the potential. Instead, I take a ^-shell potential, as in Refs. [12, US] • If the radius of this, Rq, is chosen 
to be smaller than about 0.7 fm, in the region where the waves have reached their common energy-independent forms, 
then the extracted potential will be independent of -Rn to a good approximation.* The results shown below are for 
i?o = 0.1 fm, the same radial cut-off as used in Ref. [20|, but I return to the question of the choice of cut-off at the 
end of Sec. IV, after examining the subtraction of TPE. 

To remove the dependence on the arbitrary radius Rq, I divide out the square of the asymptotic radial function, 
Eqs. © or (jlOp. from the strength of the (5-shell potential. The resulting potential is defined by 

Vs{p,r)= 1 V{p)S{r~Ro), (14) 

as in the RG analysis of Refs. [H, Equating K{p) to the DWBA matrix element of this potential, we can deduce 
the strength of the potential directly from the residual scattering amplitude, as in Ref. [20 |: 

= (15) 

IV' (p, Ra)r 

where 'ip{p, Rq) is the solution to the Schrodinger equation with the known long-range potential, and '0o(i?o) is its 
energy-independent short-distance form. 

The leading term in the short-distance interaction represents only short-distance physics, and so it should be 
independent of any low-energy (long-distance) scales such as p, or A^r. To ensure this, any dependence on these 
scales should be factored out of the normalisation of the asymptotic solutions in Eqs. Q and pH]) . In waves where 
the OPE potential can be treated perturbatively, this means normalising these functions by dividing out the energy- 
dependent factor oip^, as in Eq. ((T^). The resulting potentials in the '^F^ and channels are the defined in the same 
way as the ones used to analyse the spin-singlet channels in Ref. [20| . At LO their matrix elements are proportional 
to p^^ , showing that they are equivalent to 2L-th derivatives of 5-functions (the more conventional representations 
for the effective interactions in these partial waves). 

As discussed above, A,r should be regarded as an additional low-energy scale in the channels where OPE must be 
treated nonperturbatively. Since Plj oc l/A^, the wave functions normalised as in Eq. ([T^ still contain powers of 
A^. This can be removed by dividing out the factors of \I3lj\^^^^'^ from the zero-energy solutions. Using solutions 

with this normalisation has the effect of multiplying V{p) by A^^"''^''^. This obviously changes the magnitudes of the 
short-range interactions, but it does not affect their energy dependences. In the plots below, I just show results for 
V{p) defined using the forms given in Eq. (H^), but the extra factors would be needed if one wanted to estimate the 
momentum scales in this potential. 

Once the effects of leading OPE have been removed, the residual potential starts at order in the standard power- 
counting notation. Contributions at this order include the leading TPE potential [12, US], and a relativistic correction 
to OPE Since the power counting for the long-range forces is not affected by the renormalisation of the short- 
range ones, I denote this potential by y^^^(p). In Weinberg's counting for a partial wave with angular momentum L, 
the leading short-range term appears at the usual order, Q^^. However, if OPE is treated nonperturbatively, we need 
to take account of the fact that this term also contains a factor of A^r ■ Since A^ should be now regarded as a 

low-energy scale, we see that the net order of such a term is Q^^l"^ ^ for any L. This half-integer power agrees with 
the RG analysis in Ref. [l^, which shows that the leading term has an RG eigenvalue of 1/2. After extracting the 
effects of tensor OPE we are therefore left with an interaction that starts at order Q~^l'^ . This also contains a term 
proportional to energy (p^) at order Q^l"^ . 

Fig- m shows the short-distance interactions y'^'(p) extracted directly from various Nijmegen analyses [IJ, using 
Eq. with i?o = 0.1 fm. The partial-wave analysis, PWA93, and three potentials, NijmegenI, Nijmegenll and 
Reid93, all fit the np data with similarly good values of ^ . They can thus be regarded as alternative partial- wave 
analyses. Using results from all of them gives an indication of the systematic uncertainties associated with the different 
choices of parametrisation. This is particularly important in higher partial waves where each fit contains only a small 



* Obviously, points where the wave functions vanish in the attractive channels should be avoided for numerical reasons. Similarly very 
small values of _Ro, ~ 0.01 fm or less, should not be used in the repulsive channels. 
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FIG. 2: Plots of the short-distance interaction V^^\p), in fm^^^^, against lab kinetic energy T, in MeV. These have been 
extracted from Nijmegen PWA's or potentials using Eq.lini) with Ro = 0.1 fm, for the np channels (a) ^Pi, (b) ^L»2, (c) ^Fs, 
and (d) ^G4. 

number of parameters. In the and channels where the tensor OPE is attractive, I have taken the phase a = 
since, as already noted, their OPE phase shifts depend only weakly on a, provided the region around tt/2 is avoided. 

One immediate observation is that the residual interactions all show rapid, nonlinear dependences on energy below 
about 150 MeV. This is similar to what was found for the corresponding spin-singlet channels in Ref. |20l]- Experience 
with those channels suggests that TPE is responsible for much of this rapid energy dependence. Therefore these forces 
need to be subtracted before any conclusions can be drawn about short-range ones. 




(a) ^ (b) 



FIG. 3: Plots of the ^Po short-distance interaction V^'^\p), in fm*, against lab kinetic energy T, in MeV. Results for two choices 
of short-distance phase are shown: (a) a = 0, (b) a — 0.54. 

Finally I turn to the '^Pq wave, shown in Fig. [31 The left-hand plot was obtained using a = 0, as in the and 
cases. As in the attractive channel, this shows a rapid, non-monotonic energy dependence. In addition. 
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its overall strength is much larger than in the ^Pi channel, by a factor of 30 or more. This implies that the leading 
(order-Q^^/^) interaction in this channel is too strong to be treated in the DWBA. It should cither be iterated to 
higher orders or, possibly, treated nonperturbatively. 

It practice, it is most convenient to include this term to all orders using the equivalent short-distance parameter 
a. As discussed above, this defines a self-adjoint extension of the the OPE Hamiltonian in Eq. ([3]) by fixing the 
phase of the oscillations of the wave functions as r ^ 0. In the right-hand plot, I show the residual interaction for 
a = 0.54. This value was chosen so that the interaction vanishes at zero energy, within the uncertainties of the PWA's. 
Removing the effects of the energy-independent term in this way leaves a residual interaction that is consistent with 
a monotonic piece plus one linear in the energy. Again, more definite conclusions require subtraction of the effects of 
other long-range forces, which I now turn to. 



IV. TWO-PION EXCHANGE 



The DW method described above could also be used to separate off the scattering produced by other known 
long-range forces, most importantly those arising from two-pion exchange. However, since such forces start at order 
Q^, their effects can just be subtracted perturbatively from the residual interactions left after removal of OPE [20| . 
Although iterations of TPE will contribute at higher orders, this DWBA treatment is adequate up to order Q\ 

The relevant TPE potentials have been calculated at orders and and their forms can be found in Refs. 122. [2 3l . 
At order Q^, the long-range interactions also include a term from expanding the relativistic correction to OPE [24| . 

Vl'jir) = [y.cM + V.^{r)\ . (16) 

Lastly, there is an electromagnetic interaction, generated by 7r7 exchange, whose form is derived in [25| . 

If all of these are subtracted from K{j3) using the DWBA, then any residual long-ranged contributions to the 
scattering start at order Q^. The resulting potential is 

Note that although this potential is denoted V^^'>{p) because it contains long-range contributions of order and 
higher, it may contain short-range terms of lower order. Such terms appear first at order for P-waves in the usual 
power counting, and at order Q~^^^ for waves where tensor OPE is treated nonperturbatively. 

The long-range potentials are singular and so their matrix elements can be divergent. I therefore cut the integrals 
off, using the same radius, Rq, as in the definition of the short-range interaction. For example, the most divergent 
part of the order-Q'^ TPE potential at small radii is proportional to 1 /r^ [2^ . In a perturbative treatment of 
OPE, the wave functions behave like r-^ as r ^ 0. Together these lead to a I/Rq divergence in the P-wave matrix 
elements which can be renormalised by introducing the order-Q^ counterterm just mentioned. Other partial waves 
do not give rise to divergences to this order. 

When OPE is treated nonperturbatively the small-r forms of the wave functions are quite different, as described 
above. The r~^^^ behaviour of these leads to a divergence in the matrix element of the order-Q'^ potential, 

at least in the channels where tensor OPE is attractive. This can be renormalised by the leading short-distance 
potential, of order Q^^l"^ . In addition, there are other, weaker divergences which appear multiplied by powers of 
or m^r and which can be renormalised by higher-order counterterms. However these energy-independent terms can not 
be disentangled phenomenologically from the leading one. The relativistic correction to OPE, Eq. ([16]), also leads to 
a divergence and this can be renormalised by an order-Q'^/^ term, proportional to p^. The complete set of terms with 
orders below also includes one proportional to the square of the energy, p^. This term is of order Q^^^ and it has 
a finite coefficient at the current level of approximation. However, when long-range forces of order Q** are included, 
it will be needed to renormalise divergences from, for example, the next relativistic correction to OPE. 

The short-range terms up to order 1"^ provide a smooth, quadratic energy-dependence. Subtracting them should 
leave a residual interaction consisting only of terms of order or higher. I do this by fitting a quadratic form, 

f(^/2)(p) =Co + C2/ + C4/, (18) 

to V'^^\v) in the range T — 100 to 200 MeV. The lower end of this range is chosen to lie above the region where OPE 
and TPE can lead to rapid energy dependence. The upper limit corresponds to a relative momentum p ~ 300 MeV. 
Above this point, higher powers of the the energy could start to become noticeable, representing physics that has 
been integrated out, such as excitation of the A resonance. 
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Let me start with the two peripheral waves, where the breakdown scale for perturbation theory is so high that 
there is no question about the validity of the perturbative treatment of OPE or the standard power counting. As in 
Ref. [131, I use the full DW solutions for these, although these are not strictly necessary, to avoid the complications 
of fourth-order perturbation theory. The residual interactions for these are shown in Fig. [H In both cases, 

F^^^(p) was extracted at Rq = 0.1 fm, which is small enough that the wave functions have attained their energy- 
independent forms and that the ^F^ radial integral has converged. The integral has reached a plateau here; the 
divergences associated with the oscillatory region do not appear until i?o is less than about 0.03 fm. 

As in the corresponding spin-singlet waves [lO], the uncertainties in the PWA's make it hard to draw strong 
conclusions. The '^G4 residual interactions, and are both consistent with zero, given the spread of the 

PWA's. In the '^F^ channel, the energy dependence in V'-^^p) below 100 MeV is somewhat reduced by subtraction of 
the order-Q^'^ long-range forces. 





FIG 



4: Plots of the short-distance interaction V^'^\p), 
(b) ^G4. 



in fm^^^^, against lab kinetic energy T, in MeV, for the channels: (a) 



The interaction F'^'*) (p) in the ^Pq channel is shown in Fig.O The unsubtracted interaction, on the left, is very large, 
with a strong linear energy dependence, as a result of the divergences in the integrals of the order-Q^''^ potentials. 
The residual interaction after subtracting a quadratic fit is shown on the right."" This shows that the quadratic form 
can provide a good account of the energy dependence of ^'^'(p) in the ^Pq channel over the whole range of T up to 
300 MeV. Any contributions from forces of order and higher lie within the uncertainties of the PWA's. 




^ In this channel, I have omitted the Reid93 results from the fit since they lie systematically below the ones from the other three Nijmegen 
analyses. Including this potential would simply shift the fitted constant, and significantly increase the uncertainty associated with the 
PWA's. 
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Fig. [6] shows the results of a similar analysis of scattering in the channel. Here, the quadratic fit can account 
well for the energy dependence from T ~ 80 MeV to about 250 MeV. Below 80 MeV, the PWA's require a significant 
extra attractive interaction with a long range in view of its rapid energy dependence. It is not obvious what could be 
responsible for this, since no appropriate piece seems to be present in the order-Q* pion-cxchange potentials (26l. [27|. 
However it should be noted that its size is comparable to the uncertainties in the PWA's in this region. 



211.4 




FIG. 6: Plots of the short-distance interaction in fm", against lab kinetic energy T, in MeV, for the channel ^1)2 with 

a = 0, (a) unsubtracted, (b) quadratic fit to T = 100 — 200 MeV subtracted. 

Finally, V^^'^lp) for the ^Pi channel has a smooth, approximately linear energy dependence, as can be seen in 
Fig. [71 The value at zero energy increases as the cut-off radius Rq is decreased from 0.6 fm to about 0.2 fm. This 
is consistent with the standard power counting, where a divergence is present at order Q^. In that counting, this is 
the only divergence below order Q^. However the coefficient of the linear energy dependence in the results here also 
increases significantly over this range of Rq. This suggests that the ^Pi channel may in fact be better analysed using 
the new power counting, a conclusion that is supported by the low breakdown scale for the perturbative treatment of 
OPE found in Ref. 0. 




FIG. 7: Plots of the short-distance interaction V^'^\p), in fm^, against lab kinetic energy T, in MeV, for the channel '^Pi, (a) 
unsubtracted, (b) quadratic fit to T = 100 — 200 MeV subtracted. 

The results shown so far are all obtained with a radial cut-off i?o = 0.1 fm. This is well inside the region where 
the dependence on Rq is very small for the range of energies considered but it does correspond to momentum scales 
greater than 1 GeV, far beyond the range of validity of our EFT. This small radial cut-off was used in Ref. [l^l to 
avoid introducing artefacts proportional to positive powers of _Ro- The coefficients in the residual potential could then 
be related directly to scales of the underlying physics. The connection is more complicated here, as one would first 
need to renormalise the coefficients to remove the divergences proportional to powers of I/Rq (and there are several 
of these in the energy- independent term, involving different powers of m^r and A^r.) As a result, I do not attempt to 
make such an interpretation of the coefficients here and hence the choice of such a small radius in not crucial. 

In this case, one can ask what happens if a larger cut-off radius is chosen. The question is pertinent since Epelbaum 
and Meissner [Tsj have shown that a single extra term, as required by the new counting, can give a fairly good 
description of P- and Z?-wave phase shifts for momentum cut-offs as low as ~ 3 fm~^. Such values would also avoid 
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regions where TPE might need to be treated nonperturbatively. I have therefore repeated these analyses with larger 
cut-off radii. The results are essentially indistinguishable from the ones shown above for radii up to about 0.6 fm 
(except in regions around zeros of the wave functions). Beyond that point, the wave functions do show noticeable 
energy dependence over the range up to T = 300 MeV, as already noted from Fig. [T] This leads to a dependence on 
Rq of the residual interactions, especially at higher energies. 

If these artefacts of a finite radial cut-off are expanded as a power series in RqP^, the first three terms (up to order 
p"*) can be absorbed in the coefficients of the quadratic fit. They are therefore removed when this fitted potential is 
subtracted as described above. However, they will generate terms in the renormalised coefficients of the short-range 
potential where the momentum scale is set by 1/Ro and not by the underlying physics. This would make it difficult 
to interpret these coefficients in terms of physical scales. 

This cut-off dependence increases with energy and so it is most prominent at higher energies. A measure of how 
much the shape of the waves changes is provided by the ratio of ratios, 

^ V'(pmax,flo)-0(O,fil) ^ , . 

^ V(Pmax,i?l)V'(0,i?o) ' ^ ' 

in which the energy-dependent normalisation of the short distance wave functions cancels out. For Rq = 0.6 fm, 
Ri = 0.1 fm and Pmax = 375 MeV (corresponding to T = 300 MeV), p is greater than 0.8 for most of the waves 
considered here. The ^ 20% changes that this induces in the the residual interactions can be removed by the 
subtraction of the quadratic fit. The resulting subtracted interactions are almost indistinguishable from the ones 
shown in Figs. [5] (b) and [7] (b), except for the •^Pq which shows larger effects as a result of a nearby zero in its wave 
functions. ^ 

For Rq = 1.0 fm, the ratio p is around 0.6 for these waves. The residual interaction F'^^'(p) for the channel 
is shown in Fig. [8l The divergent terms are much smaller and so the overall magnitude of the potential is greatly 
reduced compared to the one shown in Fig. [H] (a). However the bulk of the change lies in the terms up to order p'^. 
After these have been subtracted, the differences between Figs. [6] (b) and[8](b) are small, at least for energies below 
about 250 MeV. The pattern in the ^Pi channel is similar. The ■^Pq is complicated by the fact that the wave functions 
pass through zero near 0.8 fm. Cut-off-independent results for the subtracted interaction can be obtained, but only 
for Rq ^ 0.5 fm or Rq ~ 1.2 fm. In the higher partial waves, ^F^ and ^G^, the unsubtracted residual interactions 
show very little cut-off dependence over the range 0.02 fm < Pq ^ 1 fm. 




FIG. 8: Plots of the short-distance interaction V^^\p) for the channel '^D2 with cut-off radius of 1.0 fm. Other details are the 
same as in Fig. (6] 



V. CONCLUSIONS 

Nogga, Timmermans and van Kolck [11] have found that a new power counting is needed to organise the EFT 
describing nucleon-nucleon scattering in spin-triplet channels. In particular, the leading short-distance terms in P 
and D waves are significantly promoted compared to the perturbative (Weinberg) counting. In Ref. [l^ I obtained 
a quantitative statement of this new counting by identifying the scale of the OPE potential, A^, as an additional 
low-energy scale, and then treating OPE nonperturbatively. This leads to an expansion of the short-range interaction 
describing scattering between the DW's of the OPE potential. Its terms correspond to those of an expansion of the 
DW Born amplitude in powers of the energy. The leading term is of order Q^^^^ , for any orbital angular momentum. 
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The forms of the DW's show that a nonperturbative treatment of tensor OPE, and hence the new power counting, 
is required for energies that are large enough for the waves to probe the region where the l/r^ core of the potential 
dominates over the centrifugal barrier. Otherwise the short-distance wave functions have the normal behaviour 
and perturbative counting remains valid. The analytic estimates in Ref. [l2j and the numerical wave functions in 
Sec. II both indicate that the nonperturbative approach is required in waves with L < 2 for momenta of the order of 
TOtt or larger. In F waves and above, a perturbative treatment is expected to remain valid up to energies well beyond 
the validity of the EFT. 

Here I have "deconstructed" empirical phase shifts by using DW methods to remove the effects of long-range pion- 
exchange forces. This generates a residual short-range interaction directly from the observed phase shifts. Unlike 
comparisons of phase-shift plots, this approach emphasises the low-energy region, where the EFT description ought 
to work best. The use of several Nijmegen partial wave analyses [2l| allows estimates of the uncertainties in these fits 
to the data. These can be large at low energies, implying that it may be misleading to fit the coefficients of an EFT 
to very low-energy data. 

Removing only the effects of OPE leaves residual interactions with strong energy dependences at low energies. This 
suggests that higher-order long-range interactions are also important. After removing the effects of OPE, I therefore 
use the DWB A to subtract the contributions of other long-range potentials up to order . These include TPE and 
a relativistic correction to OPE, all of which highly singular at the origin. Their DWBA matrix elements diverge as 
inverse powers of the cut-off radius but these divergences can be cancelled using counterterms at the orders required 
by the new power counting. In contrast, Weinberg counting at order would provide only one, energy-independent 
counterterm in each P wave, and none in any higher wave. 

TPE and other long-range forces up to order are able to account for much of rapid energy dependence seen 
below 100 MeV in the P and D waves. When these are subtracted the residual scattering amplitudes can, in general, 
be well fitted by three contact terms up to order Q^/^ in the new counting. The only exception is the channel, 
which seems to require an additional long-range attraction. It is not clear where this could arise, given the forms of 
the order-Q^ chiral potentials [H, HSli but one should note that the uncertainties in the PWA's for this channel are 
significant below about 70 MeV. Otherwise, TPE and the short-range forces up to order Q^^^ are able to give a good 
description of the scattering in these triplet waves up to energies of about 250 MeV. 

In the more peripheral F and G waves, the arguments of Ref. [l2| suggest that the standard Weinberg counting 
should be adequate for the energies considered here. The numerical wave functions support this, their short-range 
forms having the expected dependence driven by the centrifugal barrier. Subtraction of the long-range forces leaves 
small residual interactions, as in the corresponding singlet channels studied in Ref. [20| . Provided the cut-off radius 
is larger than about 0.1 fm, there is no sign of any divergences whose renormalisation would require the new power 
counting. 

These results indicate that the spin-triplet waves with L < 2 can be analysed consistently using the nonperturbative 
power counting developed in Refs. [HI, In addition, the results show that deconstructing scattering amphtudes, 
using the approach of Ref. [20| . can provide a very useful tool for determining effective potentials directly from 
empirical phase shifts. It should be straightforward to extend the method to coupled channels, such as '^Si-^Di, 
despite the more complicated matrix equations involved. It would also be interesting to use it to subtract two- 
and three-pion exchanges at order [2^, [23| ■ However the approach also demonstrates that there are significant 
uncertainties in the currently available Nijmegen PWA's, particularly for energies below about 80 MeV where the 
scattering is most sensitive to long-ran ge fo rces. As a result, attempts to study the importance of higher-order forces 
may require the newer PWA's of Refs. [23ll28t. when these become available. 
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